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Abstract 

We show that a fluid under strong spatially periodic confinement displays a glass transition 
within mode-coupling theory (MCT) at a much lower density than the corresponding bulk system. 
We use fluctuating hydrodynamics, with confinement imposed through a periodic potential whose 
wavelength plays an important role in our treatment. To make the calculation tractable we imple- 
ment a detailed calculation in one dimension. Although we do not expect simple \d fluids to show 
a glass transition, our results are indicative of the behaviour expected in higher dimensions. In 
a certain region of parameter space we observe a three-step relaxation reported recently in com- 
puter simulations [S.H. Krishnan, PhD thesis, Indian Institute of Science (2005); Kim et al, Eur. 
Phys. J-ST 189, 135-139 (2010)] and a glass-glass transition. We compare our results to those of 
Krakoviack, PRE 75 (2007) 031503 and Lang et al, PRL 105 (2010) 125701. 
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I. INTRODUCTION AND RESULTS 



Confined between two atomically smooth surfaces a few molecular diameters apart, fluids 
depart markedly from their bulk behaviour, displaying a prodigious increase in viscosity and 
structural relaxation times, and shear-thinning and viscoelasticity at remarkably low shear- 



adapted for shear studies 



rates and frequencies. The consensus from experiments jl-3] in a Surface Force Apparatus 



^ and simulations js-llSl is that strong confinement moves the 



system into the regime of the glass transition. Such a trend is indeed found theoretically 



14 



15| by adapting mode-coupling theory (MCT) 16l-llq| to the case of a fluid confined in a 



porous medium modelled as a random aggregate of hard spheres. MCT for a fluid between 



19| . Strictly planar walls, however. 



two smooth planar walls shows similar slowing down 
exert no force parallel to themselves, so that the dynamics remains momentum-conserving 
and long-wavelength density fluctuations travel as sound waves in directions parallel to the 



wa. 



20 



Is, unless no-slip is imposed by hand. Questions about the nature of finite-size scaling 



21| at the glass transition, and the related issue of cooperatively rearranging regions 



7| , provide further motivation for studying the effect of confinement on the slowing down of 
the dynamics of a liquid. The extraction of viscosity and friction parameters from confined 
fluid flow experiments is discussed in 



22 



23|. 



In this paper we study the dynamics of a confined dense fluid, in a coarse-grained ap- 
proach. Confinement enters the theory through an external potential; the resulting static 
inhomogeneous density background is a surrogate for the potential in our calculation. The 
mean density and temperature of the system are encoded in the static structure factor of the 



,ii9i 



m 



fluid in the absence of confinement. Our work differs in detail from those of Refs 
our use of the fluctuating hydrodynamic approach (see, e.g., encoding the interactions 
between particles and conflning medium in the free-energy functional, and examining the 
problem in detail in one space dimension. 

We summarize our main results below and in Figs. [U - [5l (i) We present a particu- 
larly transparent derivation, from the equations of fluctuating hydrodynamics, of the mode- 
coupling equations for the memory function and time-correlation function of the density 
fleld of a conflned fluid, (ii) In the fluid phase, confinement renders the density dynamics 
diffusive at long wavelengths, with a diffusivity calculable, via mode-coupling, in terms of 
properties of the inhomogeneous background density field, (iii) We show in detail (see Fig. 
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[2]) that strong enough confinement can drive the system through a glass transition, in condi- 
tions under which the system in bulk would be a fiuid. To make the calculation tractable we 
work in one dimension, where a glass transition is unlikely for a fiuid without confinement. 
However the structure of the calculation makes it clear that similar behaviour is expected 
in realistic higher- dimensional systems, (iv) The strength of potential required to produce 
the transition is lowest when the wavelength I of the potential matches the length scale 
corresponding to the structure factor peak of the fiuid (Fig. [3]). (v) For densities po below 
a value pc which depends on a continuous onset of the glassy state is observed as the 
potential strength is increased (Fig. H]). (vi) In the glassy state for po ~ we predict a 



three-step relaxation o: 



dynamics studies 



24 



the density correlation function (Fig. [S]), as seen in recent molecular 



25| of confined fiuids. (vii) Correspondingly, for a certain range of 



densities, two thresholds of confinement strength are seen for the onset of the non-ergodicity 
parameter: first a continuous onset, then a discontinuous jump. 

The prediction of a transition in one dimension, the finding that a periodic background 
suffices to enhance the transition, the emergence of diffusive dynamics for the density, and 
the multistep relaxation of the collective density correlator in the continuous transition 
regime are all features that distinguish our work from existing theoretical studies 1^ [lol of 
the glass transition under confinement. In particular, our approach offers a natural MCT- 
based explanation of the multiple plateaux in the intermediate scattering function seen 
in the simulations of {24, 25 1. The double onset of the non-ergodicity parameter seen in 
our calculations provides a much stronger example of the confinement-induced glass-glass 

n n 

transition than that seen in p^, and is similar to the behaviour found in [15[ for the tagged- 
particle density, for confinement in a disordered medium. 

The rest of this paper is organized as follows. Section [TTl begins with a discussion of 
possible dynamical effects of confinement in section III At and continues in section IIIBI with 
general conclusions from MCT for the fiuctuating hydrodynamics of confined fiuids. Section 
mil presents our calculation in detail for a model one-dimensional fluid, leading to a phase 
diagram as a function of density and confinement strength. We close in section HV] with a 
discussion and summary. 
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II. MODE-COUPLING THEORY FOR A FLUID IN A CONFINING POTENTIAL 



A. Remarks on the dynamics of confined fiuids 

Before entering into the calculations that led to the above results, some general remarks 
are necessary. The viscosity measured in the surface force apparatus (SFA) experiments 
corresponds to gradients along the z direction normal to the plates and velocity in the 

= _L plane. The density modulations contributing to these components of the viscosity 
must have a wavevector component in the ± plane, because shear with gradient along z does 
not directly affect structures varying only along z. In fact, the relevant fourier modes of 
the density must have wavevectors with nonzero z and _L components. Layering alone, with 
wavevector only along z, will not suffice. Next, what features of confinement are responsible 
for the dramatic slowing down? It seems reasonably clear that the pressure applied in the 
SFA experiments is not directly responsible, as it does not amount to a bulk hydrostatic 
compression, but simply goes into determining the thickness of the confined fluid layer. The 
confining walls rule out motion normal to their surface; they limit motion parallel to their 
surface through no-slip; and they alter the static structure of the fluid through steric or 
potential interactions. Can no-slip alone produce a significant slowing down of the density? 
To make this question concrete, consider a fluid with viscosity t] and velocity field v, with 
density field p governed by a free-energy functional F[p] (for example, a Ramakrishnan- 
Yussouff functional), lying between a pair of walls parallel to the plane, separated by 
distance w along the z direction. If we impose no-slip at the walls, the dynamics of the 
in-plane velocity field, on length scales in the _L plane large compared to w is governed 
by lubrication: {ri/w'^)\'± ~ —pV ±6F/6p. The continuity equation dtp = —V ■ (pv) then 
implies an effective two-dimensional dynamics 

S„ = V,.f^pV,^£^Uv,.fx, (1) 



T] 6p 

In writing ([1]) we have implicitly averaged over z, so that p depends only on x and y, and we 
have written the ^-averaged force density on the right-hand side in a form similar to the bulk. 



with an effective free-energy density Fpf f- Eg . ([T]) is of the type analysed in 



28 



and f is 



a multiplicative noise whose form 27|, l29|, l30| guarantees that equal-time correlations of the 
two-dimensional density field are governed by -Fg//- To linear order in density fluctuations, at 
long wavelength, ([1]) is a diffusion equation, with diffusivity kBTw^r]~^po{l — PoCq=o), where 



Cq is the direct pair correlation function of the confined fluid, and and T, respectively, are 
Boltzmann's constant and absolute temperature. Nonlinear effects enter via the interactions 
embodied in Fpf f, and can lead to a slowing down with increasing density and decreasing 



temperature 



27, 



SjJ- The confinement scale i/; in ([T]) appears in two different ways. One, 
explicitly, in determining the diffusivity ~ w'^/t], and two, implicitly, within -Fe//- The 
vanishing of the bare diffusivity oc w^, a consequence of no-slip alone, can readily be absorbed 
into a reseating of the time in ([1]), and cannot be further enhanced in an MCT feedback 
mechanism [3^. A decrease in w can lead to a glass transition only through structural 
inputs via i^e//- We are thus obliged to consider the effect of the walls on the structure of 
the fluid and thence on its dynamics [33 1. 



B. Fluctuating hydrodynamics of confined fluids 

As we remarked in the Introduction, our treatment of the effect of confinement will replace 



walls by an external potential. Consider a fluid with velocity field v and density [34] field 
p with mean po? in the presence of an externally imposed potential f/(x). The equations 
of fluctuating hydrodynamics for an isothermal fluid, extended down to the length-scales 



relevant to the slow dynamics of a fluid near structural arrest [18|, are the continuity equation 

dtp + V- (pv) = 0, (2) 
and the generalized Navier-Stokes equation 

p{dt + V ■ V)v = r]V\ + (C + r//3) VV ■ v - pV^ + f , (3) 

bp 

where C and rj are the hare bulk and shear viscosities, thermal fluctuations enter through 
the Gaussian white noise f with (f (0, 0)f (r, t)) = -2kBT[r]\S/'^ + (C + r//3) VV]5(r)(5(t). 

Let us linearise the equation of motion by ignoring 5p\ in ([2]) and ([3]), and replace velocity 
in the divergence of Eq. ([3]) using Eq. ([2]) to obtain the equation of motion for the density 
fluctuation alone 

where Dl = (C + ^v/'^)/Po- We take the space-Fourier transform of the above equation and 
obtain 



d'Sp^jt) d5p^{t) 
dt^ ^ dt 



u 



V-(pV 



6F 



zk ■ (5) 



where = Dik^ is the bare longitudinal damping coefficient. 

The confining potential U has been incorporated into the density-wave free-energy func- 
tional [26 1 

p(r) 



f3F' 



p(r) In ■ 



Po 



[p(r) - Po] 



-i / c(r - r')[p(r) - po][p(r') - Po] + /? / f/(r)p(r) 



(6) 



where = J dr, /3 = l/ksT, and c(r — r'), the direct pair correlation function in the 
absence of U, is a coarse-grained expression of the intermolecular interactions in the fiuid. 
Our reference state is m(r), the inhomogeneous equilibrium density field in the presence of 
U, which satisfies 6F^[m]/6m{r) = 0, i.e., 

m(r) 



In- 



Po 



-f3U{r) + ^ c(r - r')6m{r') 



(7) 



where 6m{r) = m(r) — po- Writing p{r,t) = po + 6m{r) + 6p{r,t), the force density from 
(E]), after replacing U{r) in terms of m(r) using Eq. (^^, takes the form 
6f3F^ 



pV 



5p(r,t) 



[5(r-r') -poc(r-r')]5p(r',t) 
V J 5m{r)c{r — r')5p{v' ,t) 

[6{r — r') — m(r)c(r — r')]5p(r',t) 



Vmfr) 



m r 



-Sp{r,t)V c{r-r')5p{v',t). 



(8) 



Note that the effect of confinement is contained only in the background density field m(r); 
U itself does not appear explicitly. The hydrodynamic equations and ([3]) with (IHl) then 



readily yield the dynamical equation 



35| 



-tk ■ (Fk + fk) = ■ (Fr + Fr + fk^) 



(9) 



for the spatial Fourier transform 5pk(t) of the density field, with force densities 

F^%t) = tksT I [kcq + (k - q)/po^i°)]5mk-q5pq(t) (10) 

arising from interaction with the static inhomogeneous b ackg round, to first order in (5m(r), 
which suffices for the one- loop treatment we will present 36|, and 

F^''(t) = '-ksT j [qcq + (k - q)ck_J5pq(t)5pk-q(t) (11) 



from the pairwise interaction of fluid density fluctuations. In (Q fk (^) is the longitudinal 
part of the Fourier-transform of the bare noise f in ([3]), and S^'^ in iQ and ffTOl) is the static 
structure factor of the bulk fluid without confinement. In (ITU]) and fllip and hereafter 
denotes J d'^q/{27iY. The autocorrelation of f^{t) is linked to the bare longitudinal damping 
Fk = Di^k^, and the Kubo formula [s^] tells us that the excess damping due to interactions, 
expressed in the time domain, is given by the memory function 

A^k(t) = ^^(Fk(0)-F_k(t)) (12) 

where ^ — )■ oo is the system volume, and Fk(t) is as in (p|)- (fTT|) . 

We define translation-invariant correlation functions S'k(t) = / (i"'r exp(—zk • 
r)[(5p(ro)(0)5p(ro + r)(t))]r(, where [ ]ro denotes an average over tq over a period of the 
background. Within a Gaussian decoupling approximation (l9i)-( |T2l) lead to 



k{t) + DLk^^{t) + ^^M) + / Mk{t - T)MT)dT = (13) 

Jo 

for Mt) = ^k(t)/5k(0), with 

Mkit) = ^ y"[k ■ qc, + k ■ (k - q)ck^,]^Sk^^{t)S^{t) 

+^ ^[k . qc, + k ■ (k - q)/po?Sl_^S^{t), (14) 

where Ci = ksTpQ and S*^ = J (i''rexp(— zk ■ r)[5m(ro)5m(ro + r)]ro is the structure factor 
of the background density field. Eqs. (fT3l) and ((Hj) are the equations of mode-coupling 
theory (MOT) for a confined fluid, as an initial value problem for S'k(t) given the static 
structure factor ^k = 5'k(t = 0) of the confined fluid. These equations are identical to those 
obtained by Krakoviack [l^ through the projection operator formalism for the case of a fluid 
confined by a porous medium. Our approach offers a simple and transparent derivation of 
these results. 

Irrespective of the details of the confining medium, ( 1141) leads us to our first result, an 
important general feature of confined MCT. Of the two contributions to A^k(^) on the right- 
hand side of (fT4l) . the second, coming from the interaction of dynamic density fluctuations 
with the background density field, is non-vanishing for wavevector k — )■ 0. The interac- 
tion of the fluid with inhomogeneities of the confining medium damps the flow even at 
zero wavenumber, and MCT provides a convenient way to calculate this effective "Darcy" 
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damping. It is this contribution that leads to an effective no-shp condition on confining 
walls endowed with structure, periodic or otherwise. As a consequence, for long timescales 
and small k the third and fourth terms on the left-hand side of fll3p dominate, reducing 



it, as is to be expected 



39, 



40|, to a diffusion equation for 0k (t) with collective diffusiv- 



ity ksT/ S]f_=o-M. 00 where A^oo is the Fourier-transform of the memory function at zero 
wavenumber and frequency. The crossover to Brownian rather than Newtonian dynamics 
assumed for convenience in fl^. l4l| 42 1 emerges at small wavenumber in a calculable form 
as a result of interaction with the confining medium. 



III. MODE-COUPLING IN ONE DIMENSION IN A PERIODIC POTENTIAL 

The theory presented in the previous section is quite general and can, in principle, be 
applied in any dimension. The mode-coupling calculation for a system with spatially periodic 
confinement in two or three dimensions makes enormous demands on computational time, 
because of the loss of isotropy. We therefore chose to display the principle of the calculation 
by working in one dimension. We are aware that a simple fluid of, say, hard rods in one space 
dimension has no crystalline phase and therefore cannot be supercooled or overcompressed, 
rendering it a poor candidate for a glass transition. As MCT is relatively insensitive to 
spatial dimensionality, we expect nonetheless that the trends in our Id MCT treatment will 
be found in calculations in dimensions > 2. We will return to the question of the behaviour 
of one- dimensional systems at the end of the paper. 

In order to implement our calculation we need an expression for for the confined 
fluid. Rather than taking it from experiment or liquid theory, we express it in terms of the 
structure factor S^^ of the bulk fluid which we treat as given. We expand (|6]) around the 
static inhomogeneous density m(r); S'k is then determined by the coefficient of the term 
quadratic in 6p. This crude approximation is adequate for the purpose of illustration, and 
can be improved upon if necessary. To linear order in (5mk it is straightforward to see that 

5k = 4°^ + ^mk^f . (15) 

Po 

We must solve the closed set of equations f|T3l) - f|T5l) numerically to obtain the behaviour 
of our confined fluid. As remarked above, we do this for space dimension d = 1. The 
calculation can be viewed as a schematic MCT for the confined fluid problem in which the 
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FIG. 1: (Color online) (a) A schematic presentation of a one-dimensional fluid confined by periodic 
external walls, which we re-express in terms of a periodic mean density profile, (b) The direct 
correlation function and the static structure factor 5^*^^ for a one-dimensional bulk fluid, (c) 
The Fourier transformed background density Jm^. (d) The bulk structure factor 5^*^"* and its 
modification Sk due to the external potential. 

structure of the confining walls is taken into account in a simple manner. 

We first solved the MCT equations (fT3|) and (fT4|) without confinement, in one dimension, 



using the direct correlation function for a one-dimensional hard- rod fluid [43| as input, and 
found an MCT glass transition at po = 0.7726 (Fig. [2] (a)). The nature of the MCT 
calculation in d = 1 deserves some explanation. Through the quadratic nonlinearity, one 
mode at wavevector k couples to two others at ki, \^2 with ki + k2 = k. For c/ > 2 all three 
of these modes can be taken to lie on the first shell of maxima |k| = of the structure 
factor ^(k). How then do we understand the MCT glass transition qualitatively for d = 1, 
where the "shell" is two points? The answer: for a Id fluid the higher-order peaks of 5'(k) 
carry substantial weight; the dominant triple of modes coupled by MCT must consist of two 
with I A; I = fco and one with \k\ = 2ko. Having established an MCT glass transition in d = 1, 
we now proceed to examine its modification by confinement. 

Consider one-dimensional confinement by a periodic potential U{x) = U sin^, in units 
of /cfiT, where the parameters U and i control the strength and periodicity respectively of 
the potential. Given u{x) we can construct the background density field 6m[x) = m{x) — po 
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FIG. 2: (Color online) (a) The relaxation of the normalised coherent intermediate scattering func- 
tion <^k(i) for k = 5.0 for an unconfined one-dimensional fluid shows the mode coupling transition 
at density po = 0.7726. (b) We set the density po = 0.70 so that the system is in a fluid state 
[(j)]i{t) relaxes to zero] in the absence of an external potential. For an external potential with period 
i = 1.3, corresponding to first maximum of the structure factor of the fluid without confinement, 
we increase M from a small value and observe the MCT transition at M = 12.2 x 10~^. The plot 
shows the relaxation of (j)]ii't) for k = 5.0. 

as in Eq. ([7]). Instead of specifying u{x), we will therefore characterize our confining medium 
by specifying 6m{x), which we take for simplicity to be sinusoidal, 6m{x) ~ Msin27rx/^. 
We set the density po = 0.70, for which the system in the absence of a confining potential is 
in the fluid state, and take i = 1.3, which means the period of the external potential is 27r 
over the wavenumber at which the structure factor of the bulk one-dimensional fluid has its 
primary peak. As M is increased past a threshold of about 12.2 x 10"^ we find, solving ( !T3|l 
and f[T^ with f[T^ . that the mean relaxation time obtained from 0k(^) diverges, and the 
non-ergodicity parameter = \imt-^oo 0k(^) jumps to a nonzero value (Fig. [2] (b)). This is 
the Id MCT glass transition induced by confinement, i.e., by a periodic potential. Keeping 
Po fixed, we ask how changing £ affects the threshold M for the transition. In Fig. [3] (b), we 
scan around the first peak of the structure factor and find that at high density, the threshold 
value of M is the lowest, i.e., the strongest enhancement of the transition is achieved, at a 
value of £ = 1.3 corresponding to the peak. Values of £ on either side of 1.3 require a larger 
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threshold value of M. However, at low density, the threshold of M is very insensitive to 
I. This can be understood as follows. At higher densities, within MCT, the particles form 
a cage and i = 1.3, being compatible with the preferred interparticle distance, facilitates 
this caging. For other values of i the two length scales are different and, hence, 6mk is less 
effective. At low density, there is no caging, so varying i has little effect. It is important to 
note that the transition at po = 0.60 is always continuous whereas at po = 0.70 it can be 
continuous or discontinuous depending on the value of i (Fig. [3] (d)). 

We list several features of our phase diagram that are consistent with the findings of 
[l^ for fluids in disordered porous media. For moderately high density po, the effect of the 
external potential can be seen as enhancing the correlations already present in the fluid, 
i.e., promoting the effect of the quadratic term in the MCT equation f[T^ . At substantially 
lower densities, where correlations in the bulk fluid are weak, it is still possible within MCT 
to get a transition to a non-ergodic state by increasing M, but now the dominant role is 
presumably played by the term linear in Sq(t) in ( fT4l) . Correspondingly, for densities po lower 
than a value pc which depends on i, the transition turns continuous (Fig. HI). Quantitatively, 



however, the medium has a weaker effect at high densities than in iJ]. More puzzling, it 
is substantially more effective at low densities than in UM. Re-entrance is another feature 
observed in common with pLj]: at densities somewhat below pc the threshold value of M 
decreases with decreasing density. It is surprising that we observe it in our simple one- 
dimensional periodic system, and we suggest a simple interpretation. At low densities, 
the glassy state, if it exists, presumably comes from a feedback-enhanced slowing down of 
single-particle crossings of barriers posed by the external potential. As density is increased, 
interparticle repulsion becomes more important, and particles typically occupy regions not 
near the minima of the external potential. Thus, repulsion lowers the effective kinetic barrier 
to motion. Possibly our explanation is related to that proposed in |l4|, but this is unclear. 

On the continuous MCT-glass transition line, there are possible difficulties with infrared 
divergences in the mode-coupling integral ([HI) because we are working in one dimension. 
This issue arises in a pronounced manner for a single particle in a disordered medium 44| . 
where there is no transition in d = 1 because the particle is always localized. There are 
however important differences between our system and that of 4^. Our model is probably 
free from the /c — )■ problem because the background medium in our case is periodic, not 
disordered, and lacks weight at k = 0. We have nevertheless checked that varying the lowest 
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FIG. 3: (Color online) (a) The static structure factor Sj^ for an unconfined one-dimensional fluid 
at a density pQ = 0.70. The various I values marked on the fc-axis in the figure correspond to the 
peak position of 5mk [Fig. [Il^c)]. (b) The phase diagram for M vs £ for two different densities. 
The threshold of M is quite sensitive to i for pQ = 0.70 whereas it is less sensitive at po = 0.60. (c) 
Phase diagram for po vs M for four values of i. (d) As in (c), for a larger range of pQ. In both (c) 
and (d), filled and open symbols respectively indicate continuous and discontinuous transitions. 
Note: For clarity, the extension of the discontinuous transition line beyond the point where it 
crosses the continuous transition is not shown here [see Fig. [5] 

(a)]. 
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FIG. 4: (color online) At density po = 0.60 and i = 0.80, the non-ergodicity parameter = 
limf_j.oo 4'k.{t) shows a continuous onset to the MCT glass transition with increasing strength of the 
external potential. The plot shows the relaxation of <j)k{t) for k = 5.0. 

wavenumber over a limited range (from 0.1 to 0.2) does not alter, to 0.1 % accuracy, the 
MCT transition density as estimated by the onset of a time-persistent density correlator. 

For po moderately large, but still below p^, fk undergoes two onsets as a function of 
M. Across a first threshold value of M, fk rises continuously from zero, presumably driven 
by the linear term in the MCT equation. Upon increasing M past a second threshold, a 
disconfnuous jump ■„ /. is seen, which amounts to a glass-glass Uans.fon |F,g. 0a)l. A 
hint of such behaviour has been reported in [IJ] ; the effect appears far stronger in our case. 
The relaxation of the normalized intermediate scattering function in this second glassy 
state [see Fig. [5](b)] shows two plateaux followed by a third nonzero asymptotic value. 
Such dynamical behaviour has been reported for the self-intermediate scattering function in 
molecular dynamics studies of Lennard- Jones fluids under planar confinement to a thickness 



of about three molecules by structured walls 2^], where it was called three-step relaxation. 



More recently similar dynamics has been reported for the coherent part of the intermediate 



scattering function for hard spheres in random media 25|]. The possibility of such multistep 



relaxation scenario was predicted in the context of various possible MCT integrals 



45 



and obtained by Krakoviack 15|] for tagged-particle motion, from MCT in a disordered 



medium. Such a complex relaxation scenario has also been reported in a variety of other 



contexts 



47^49] 
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FIG. 5: (Color online) (a) For 0.65 < po < 0.72, i = 0.8, a continuous transition from liquid to glass 
(segment BD) is followed by a discontinuous glass-glass transition (segment BC). Note, however, 
that we do not resolve the phase boundary close to point B. (b) If the density pQ is close to but lower 
than />c, the threshold density where the transition switches from discontinuous to continuous, the 
density correlation function shows a three-step relaxation for strong enough background density. 
The figure shows the relaxation of (p^^^it) at a density pQ = 0.70 as a function of log(time) as we 
increase the strength of background density with I = 0.80. 

IV. SUMMARY AND DISCUSSION 

We have succeeded in providing an economical description of the slowing down of a fluid 
under confinement. We derived the mode-coupling equations through the dynamic density- 



wave approach [50| of fluctuating hydrodynamics, and implemented confinement in the form 
of an externally imposed periodic potential. In the spirit of a one-loop approach, we retained 
interactions between the fluid and the background density only to bilinear order, but this is 
not a serious difficulty; we could work with the background molecular field instead. In order 
to make our numerical calculations manageable we chose to work with a one-dimensional 
model and to deal with properties averaged over one period of the potential. Once we relax 
the period-averaging constraint, we can predict the degree of slowing down as a function 
of location in the potential. This is the nearest analogue, in our Id model, to calculating 
properties as a function of distance from confining walls. Remarkably, our model calculation 



14 



reproduces all the features observed in more detailed treatments, including a crossover to a 
continuous transition, the phenomenon of re-entrance, at lower densities, and a glass-glass 
transition as a function of confinement strength at intermediate densities. Unlike in 14| . 



quenched disorder plays no role, and the density dynamics is diffusive at small wavenumber, 
in contrast to 19|. In future work we will consider in detail the problem of confinement in a 



planar geometry with structured walls, thereby improving on the treatments presented here. 

Apart from the glass-glass transition, which appears as a prominent feature in our treat- 
ment unlike in 1^ where it is barely detectable, our most novel predictions are (i) that even 
a one-dimensional fluid, at least within an MCT approach, can undergo a glass transition; 
(ii) that an imposed periodic potential causes this transition to occur at lower densities and 
higher temperatures; (iii) that in a calculable range of densities, the relaxation of the density 
correlation function should take place through a three-step process, an effect for which there 
is evidence in molecular dynamics studies of two types of confined fluid systems 



and as shown in MCT with confinement in a disordered medium [151] for the tagged-density 
correlation function. 

The initial experimental observations of the effect of confinement were done in shear 
flow measurements. It is therefore imperative to extend our approach to include shear. We 
lave done this in a simplified implementation of the ideas of 5l|. Our preliminary results 



52| on an isotropized version of the calculation find that the shear-thinning moves to very 



low imposed flow-rates in the presence of the confining potential. Further results require 
imposing shear with planar confinement, preferably with structured walls. We are confident 
that our results will change only quantitatively when these elements of greater realism are 
introduced. 

Lastly, should one expect a glass transition, or at least an MCT-style enhancement of 
viscosity, in experiments or simulations on one-dimensional fluids, confined or otherwise? 
It might seem unlikely, given that a, Id system has no crystalline phase past which one 
can supercool it to enter a region of metastability and glas siness. However, supercooling is 
not necessary for viscosity increase of the MCT variety 53|, and numerical calorimetry [5^ 
finds precursors of glassy behaviour even for Id Lennard- Jones mixtures. It might therefore 
be worth testing our ideas in simulations on Id fluids, perhaps multicomponent or with 
complicated interactions, with and without confinement. In any case, we look forward to 
tests of our predictions in experiments and simulations, whether in one or higher dimensions. 
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